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ABSTRACT 

We report systematic variations in the emission line ratio of the CO J = 2 — 1 and J = 1 — 
transitions (i?2-i/i-o) in the grand-design spiral galaxy M51. The i?2-i/i-o ratio shows clear evidence 
for the evolution of molecular gas from the upstream interarm regions, passage into the spiral arms 
and back into the downstream interarm regions. In the interarm regions, i?2-i/i-o is typically < 0.7 
(and often 0.4-0.6); this is similar to the ratios observed in Galactic giant molecular clouds (GMCs) 
with low far infrared luminosities. However, the ratio rises to > 0.7 (often 0.8-1.0) in the spiral arms, 
particularly at the leading (downstream) edge of the molecular arms. These trends are similar to those 
seen in Galactic GMCs with OB star formation (presumably in the Galactic spiral arms). i?2-i/i-o 
is also high, ^ 0.8 — 1.0, in the central region of M51. Analysis of the molecular excitation using a 
Large Velocity Gradient radiative transfer calculation provides insight into the changes in the physical 
conditions of molecular gas between the arm and interarm regions: cold and low density gas (< 10 
K, < 300 cm~^) is required for the interarm GMCs but this gas must become warmer and/or denser 
in the more active star forming spiral arms. The ratio i?2-i/i-o is higher in areas of high 24/im 
dust surface brightness (which is an approximate tracer of star formation rate surface density) and 
high CO(l-O) integrated intensity (i.e., a well-calibrated tracer of total molecular gas surface density). 
The systematic enhancement of the C0(2-l) line relative to CO(l-O) in luminous star forming regions 
suggests that some caution is needed when using C0(2-l) as a tracer of bulk molecular gas mass, 
especially when galactic structures are resolved. 

Subject headings: galaxies: individual (NGC 5194, M51) - ISM: clouds - ISM: evolution 



1. INTRODUCTION 

The standard, albeit simplistic, picture of interstellar 
matter (ISM) phases posits that giant molecular clouds 
(GMCs) are assembled in spiral arm shocks from diffuse 
interarm HI gas and then photo-dissociated back into 
the atomic phase by OB star formation within the spi- 
ral arms. This picture predicts a rapid gas-phase change 
across spiral arms from atomic to molecular and back 
into the atomic after spiral arm passage. New CO(l-O) 
observations of M51 suggest a very different picture in 
which there is little (or no) gas-phase change between in- 
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terarm and arm regions with the majority, 70-80%, of the 
neutral gas remaining molecular from arm entry through 
the inter-arm regio n and into the next spiral arm passage 
(jKoda et al.ll2009[ ). The molecular gas fraction could be 
even highe r if dark H2 gas exis ts (e.g., Planck Collabora- 
tion 2011; iGrenier et al.|[2005l ). which is not detectable in 
CO emission. In this new picture, the obvious molecular 
spiral arms are assembled from pre-existing small GMCs 
by spiral arm forcing and orbit crowding within the arms. 
A fraction of the gas forms dense cores i n the molec- 
ular spiral arms, leading to star formation (|Egusa et al.l 
l2011f ). but the majority is shredded apart by spiral shear- 
ing motions an d goes back to the intera rm regions as 
smaller GMCs (IKoda et al.l [20091: 1 Wada fc Koda 2004) . 

Questions still remain as to what the physical condi- 
tions of the abundant molecular gas in the interarm re- 
gions are and how these conditions change within the 
spiral arms. The lowest rotational transition J = 1 — 
of the carbon monoxide (CO) molecule is the most reli- 
able tracer of the overall molecular gas mass and surface 
density in galaxies. The transition energy of the CO(l-O) 
corresponds to E/k ~ 5.5 K and this is conveniently be- 
low the typical observed temperatures of resolved GMCs 
(^ 10 K; Scoville & Sanders 1987); in addition, the crit- 
ical density for CO(l-O) collisional excitation in the pres- 
ence of photon trapping s a few x 100 H2 cm~"^ which is 
roughly the same as th e average density within the GMCs 
(|Solomon et al.l[T987l ). Therefore, CO(l-O) line emission 
arises from the levels which are close to thermalized, 
and the C O (1-0) luminosity traces their bulk molecular 
mass fsee lScoville fc Sanders! [l987h . Empirically, C0(1- 
0) luminosities of GMCs are correlated with their virial 
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mass estimates. A nearly constant mass-to-light ratio 
has been calibr ated with virial masses of resolved GMCs 
in the Galaxy ([Solomon et al.l Il987t 'Scoville & Sanders] 
19 87') and in nearby spiral g alaxies |Bolatto et al. 200a 
Donovan Mever et al.l 12013 ). [Note that, despite the 
constant ratio observed in normal spiral galaxies, there 
are discu ssions on its variations in l ow-me tallicity dwarf 
galaxies ([Israeli fl99% iLerov et ahl l2009b[) . where self- 
shielding of GMCs from ambient ultraviolet radiation 
may not be as effective.] 

Within the GMCs only a small fraction of the gas 
is typically within the dense core regions leading to 
star formation. HCN(l-O), having a higher critical den- 
sity (- lO^cm-^; iKohno et ahl 119961 ) than low- J CO 
lines, traces the emission from such dense cores and 
this line emission exhibits a linear correlation with in- 
frared (IR) emission, a tracer of embedded star forma- 
tion, over 7-8 orders of magnitude - from starburst galax- 
ies down to small dense cores within Galactic GMCs 
([Gao fc Solomonll200llWu et al.ll2?)05h . Recentlv, C0(3- 
2) ha s also been used to tr ace dense cores (iMuraok a et al.l 
120071: iWilson et all |2009[) . Its higher critical density 
(> lO^'cm"^) and relatively high temperature corre- 
sponding to the upper energy level (~ 33 K) imply that 
it traces the environments with active star formation - 
either dense star-forming molecular cores or gas heated 
by recent star formatio n. This has been demonstrated 
by iKomugi et al.l ([20071 ) who found a linear correlation 
between surface densities of CO(3-2) emission and star 
formation in nearby galaxies. 

C0(2-l) emission has an upper-level energy temper- 
ature of ~ 16.5 K and a collisional critical density 
of ~ 10^~'*cm~'^, slightly above the standard CMC 
temperature and density; hence the CO(2-1)/CO(1-0) 
line ratio (-R2-1/1-0) can trace the change of phys- 
ical conditions of the bulk molecular gas. In fact, 
the ratio shows a systematic variation in the Galaxy 
([Sakamoto et al.l [199 9) and has been used to investi- 
gate the varying physic al conditions ([Sakamoto et al.l 
Il997t ISawada et al1l2001l ). i?2-i/i-o has also been used 
to infer gas physical conditions in external ga l axies , 
especially in their central regions (IKnapp et al.l [19801: 
Braine fc Combes! Il992i: iTurner et al.lll993l: lAalto et al.l 
19951 : iHarrison et al. 

Attempts to map j?2 -i/i-o across galaxy disks have 
remained inconclusive. IGarcia-Burillo et al.l ([1993D an- 
alyzed their 2-1 and 1-0 imaging of M51 and found no 
systematic variation in _R2-i/i-o between the spiral arm 
and interarm regions. However, the early CO(l-O) data 
([Nakai et al.lll994) used in their analysis appear to be in- 
consistent by a factor of ^ 2 in flu x with four other later 
measurements ([Koda et al.ll2011[). perhaps due to cali - 
bration difficulties 20 years ago. ISchinnerer "eFall ([2OIOI) 
analyzed i?2~i/i-0i as well as other molecular lines, in 
the western spiral arm of M51 using data obtained with 
interfer ometers, which was unfortunately combined with 
the old iNakai et al.l ([1994f ) single-dish data with the ap- 
parent calibration problem. Therefore, the question of 
large-scale variation of i?2_i/i_o remains unanswered. 
Here, we present new C0(2-l) and CO(l-O) imaging with 
a new and consistent observation technique and care- 
ful calibration for both lines. These new data permit 
much-improved measurements of the spatial variations 



in i?2-i/i-o in both arm and interarm regions of M51. 

2. CO AND INFRARED DATA 

The CO J = 1 — data are taken from iKoda et al.l 
(Unni), which were taken as a part of the CARMA 
NObeyama Nearby-galaxies (CANON) survey. The orig- 
inal observations include both interferometer data from 
the Combined Array for Research in Millimeter Astron- 
omy (CARMA) and single-dish data from the Nobeyama 
Radio Observatory 45m telescope (NR045). However, 
we use only the NR045 data here to avoid unnecessary 
concerns, if any, about potential imaging artifacts intro- 
duced in the process of image reconstruction. CARMA 
data, observing small angular scales, do not play a role 
in our analysis, since we use a spatial resolution of 19.7". 

The improved accuracy of the new data is due 
to the 25-BEam Array Receiver System (BEARS; 
'Sunada et al.' ^OOOl and the On-The-Fly (OTF) mapping 
technique (Sawada et al. 200j). They enable us to map 
the large galaxy within a short observing duration, en- 
suring consistency in calibration. Multiple OTF scans in 
the RA and DEC directions are averaged to minimize sys- 
tematic errors which appear dominantly along the scan 
directions (note that this type of error could not be char- 
acterized without the OTF technique). Therefore, the 
new data provide a much higher accuracy in relative cal- 
ibration across the map. The absolute flux calibration 
is consistent with other independent measurements at a 
few perc ent level. T h e det ails of the data reduction are 
given in iKoda et al.l ([20111 ) . The final data cube has a 
spatial resolution of 19.7" (degraded by gridding from 
the beam size of 15") and a Icr noise of '--^ 37 mK on the 
main beam temperature scale in a lOkms"^ channel at 
a pixel size of 5.96" x 5.96". It covers the entire molec- 
ular disks of NGC 5194 (main galaxy) and NGC 5195 
(companion). 

The CO J = 2 - 1 data are from lSchuster et al.l([2007D . 
which are distributed as a part of t he HERA CO Lin e 
Extragalactic Survey (HERACLES; ILerov et al.l[2009al ) . 
The data are taken with the 18 element focal plane het- 
erodyne receiver array (HERA) on the IRAM-30m tele- 
scope. The OTF technique was also employed. The 2-1 
data cover the entire disk of NGC 5194, but do not in- 
clud e NGC 5195. The orig inal integrated intensity image 
is in lSchuster et al.l ([20071 ). The spatial resolution is 11", 
and the noise is ~ 22 mK in a lOkms^^ channel on the 
main beam temperature scale. The C0(2-l) data cube 
was smoothed to the 19.7" resolution and regridded to 
the CO (1-0) reference grid. When comparing the two 
maps, we noticed that the global C0(2-l) distribution is 
offset by ~ 1 pixel in the south-west direction with re- 
spect to CO(l-O). We therefore shifted the C0(2-l) data 
correspondingly before making the integrated intensity 
maps with both lines integrated over the same velocity 
ranges referenced to the observed rotational velocity field 
in the galaxy. This treatment does not change the con- 
clusions in this paper, as the spiral arms and interarm 
regions are considerably wider than one pixel. 

The Spitzer 24/ito image and Ha image are taken from 
the archive of the Spitzer Infrared N earby Galaxies Sur- 
vey (SINGS; IKennicutt et"all [2 00l). The GALEX far- 
uhraviolet (FUV; 1539A) image ([Gil de Paz et aLllMl 
is taken through GalexView. We adopt the central co- 
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Fig. 1. — Maps of the central 5.4' X 4.1' region of M51 in the J2000 coordinate system, (a) CO J = 1 — integrated intensity, /cO(i— 0)) 
map in units of main beam temperature, (b) CO J = 2 — 1 integrated intensity, /co(2-l)i map in units of main beam temperature, (c) 
CO J = 2 — 1/1 — Hne ratio, -R2— i/i— o (= ^CO(2— l)/^CO(l-0) )i map. (d) Spitzer 24/im surface brightness, /24^m> image. The signal to 
noise ratios are high across the entire area including at the edges, i.e., > lOcr in both CO(l-O) and CO(2-l). The circle at the bottom-right 
corner in (a) is the beam size of the CO J = 1 — map, 19". 7 {~ 780 pc at the distance of 8.2 Mpc). All other images are smoothed to 
this resolution. Contours in all panels are at /cO(i— 0) = 12, 24, 36, and 48K • kms~^. The line ratio varies systematically between the 
spiral arms and intcrarm regions. 
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ordi nates of (a, J)j2ooo = (1 3h29m52.7s, +47dllni42.8s) 
from iHagiwara et all (|200lD . 



3. CO LINE RATIOS 

Line ratio images require high significance detections, 
especiahy of the denominator emission hne; therefore, 
in the outer disk of NGC5194 where the hne bright- 
ness decreases, smoothing and/or chpping of low signif- 
icance data is required. We adaptively smoothed the 
C O (1-0) data with the s patial binning method developed 
by iCappellari fc Copinl ((2003) and used data clipping in 
these areas when a large number of pixels (> 10) need 
to be averaged; identical smoothing and clipping was ap- 
plied to the CO (2-1) map. Figure [TJi and b show the 
resuhant maps for CO(l-O) and C0(2-l). 

The noise levels change across the maps as a result of 
the variable velocity widths for integration and adaptive 
smoothing, but the signal to noise ratios are high across 
the entire detected area, including at the edges, i.e., > 
lOcr in both CO(l-O) and C0(2-l). 

The 24/^m image is processed in the same way and is 
shown in Figure [TJl. The 24/xm emission appears at the 
leading (downstream) edges/sides of the CO(l-O) spiral 
arms (contours). Figure show the Ha and FUV 
images at higher resolutions. These emissions, tracing 
recent star formation, are in general at the leading edges 
of the molecular spiral arms. 

3.1. Spatial Distribution 

Figure [1]: is the -R2-1/1-0 niap and shows a spiral pat- 
tern of elevated i?2-i/i-0j i-C-j -R2-1/1-0 varies system- 
atically between the spiral arm and interarm regions. In 
the interarm regions, this ratio is < 0.7 and often as 
small as ^ 0.4 — 0.6, while around the molecular spi- 
ral arms it is elevated to > 0.7 and often becomes as 
high as ^ 0.8 — 1.0, mainly at the leading (downstream) 
edges of the arms (see also Figure [2^, b) . The average 
i?2- 1/1-0 over the entire disk is ~ 0.70. High ratios also 
appear occasionally at some upstream edges (e.g., in the 
north-west quadrant). We note again that both CO(l-O) 
and CO (2-1) are detected significantly even at the edge 
of this map. In the central 70" (^ 2.8 kpc), i?2-i/i-o is 
also high, ~ 0.8 — 1.0. The high ratio appears to coincide 
with strong 24^m emission (Figure [TJi). Figure [3] shows 
the phase diagrams of the CO(l-O) integrated intensity 
and i?2-i/i-0i again evidencing the association of the 
high ratio with the molecular spiral arms, especially at 
their leading edges/sides (i.e., right-hand side). 

The spatial variations of i?2-i/i-o across spiral arms 
are also confirmed in Figure |4l We adopt a simple log- 
arithmic spiral pattern with a pitch angle of 20 deg and 
measure statistical quantities around the spiral curves 
over 40 deg segments (20 deg increment). The average 
i?2- 1/1-0 increases to > 0.7 in the molecular spiral arms 
(blue), especially on/near their leading edges (red), and 
decreases to < 0.7 in the interarm regions. The areal 
fraction of high ratio gas (0.8-1.0) is high in the spiral 
arms and their leading edges, but low in the interarm 
regions. The fraction of low ratio gas (0.4-0.6) is low in 
the arms, but high in the interarm regions. The central 
region also shows a high ratio. The average is 0.78 within 
the central 90" (-- 3.6 kpc) region and 0.82 within 70" 
(~ 2.8 kpc). [Example histograms of i?2-i/i-o in the 



regions are presented in Appendix [X] and confirm the 
significant difference among the regions.] 

Figure [2^ shows the distribution of gas with low 
(i?2-i/i-o = 0.4-0.6), medium (0.6-0.8), and high ra- 
tios (0.8-1.0) with respect to the molecular spiral arms 
(contours). Figure [2)d highlights the distribution of the 
high ratios (0.8-1.0). For reference, this figure includes 
the outer areas where a large number of pixels (> 10) 
need to be averaged. The variations of i?2-i/i-o on large 
scales are evident: high ratios around the molecular spi- 
ral arms (contours) as well as in the center, and lower 
ratios in the interarm regions. The high ratios appear 
mostly at the leading edge/side of the molecular spiral 
arms where OB stars are predominantly located (Figure 

m, d). 

3.2. Correlations 

Figure [5] shows pixel-by-pixel correlations of i?2-i/i-o 
with three parameters (for the pixels in Figure [1]). The 
24^m surface brightness /24pm is roughly a proportional 
of the star formation rate surface density (jCalzetti et al.l 
I2007D and the CO(l-O) integrated intensity /co(i-o) is 
proportional to the molecular gas surface density. Their 
ratio f24:fim / Icoii-o) is a measure of the star formation 
efficiency. [Note that ICalzetti et al.l(j2007D discussed that 
the combination of 24/xm and Ha emission provides a 
more accurate star formation rate; however, f2ip,m is 
used here as an approximate tracer, since it shows in- 
herent correlations. Correlations with some other star 
formation rate tracers are given in Appendix IbI] 

Most remarkable are the lower envelopes, or cutoffs, of 
the distributions in Figure[5](left panels). No data points 
exist in the lower-right corners. We note that there are 
some dependences between x and y axes in the middle 
and bottom panels (though not in the top panels) as both 
axes include /co(i-o)- However, we conclude that the 
cutoffs in all panels are intrinsic to physical conditions 
of the gas, since our detection limits cover the entire 
areas in these plots. Hence, _R2-i/i-o increases with the 
intensity of star formation and gas surface density. 

The right panels show correlations using averages and 
their scatter. -R2- 1/1-0 increases clearly with star for- 
mation activities (i.e., /24/im and f24fim/ Ico{i~o}) and 
weakly with the gas surface density (/co(i-o))- The 
weaker correlation is due to the spatial offset between the 
molecular spiral arm and high i?2-i/i-o ( ^3.ip . Multiple 
physical processes may contribute to the intrinsic scat- 
ters in the left panels, and therefore, linear model fits 



may be inappropriate. For instance, R 



2-1/1-0 



can be 



enhanced due to increases in density prior to star for- 
mation and increases in gas temperature due to heating 
after star formation (Q. 

In addition to the clear correlations with star forma- 
tion activities, some high values of i?2-i/i-o are seen 
with low /24Mm: ^co(i-o), and /24A«n/^cO(i-o)- FigurelH 
shows the locations of this component, i.e., with a high 
ratio, -R2-1/1-0 > 0.8 and low star formation efficiency, 
/24Mm/^co(i-o) < 0.6, with respect to f24p.m (white con- 
tours) and Ha emission (right panel) . Some regions with 
the elevated ratios appear at the upstream edges of the 
24/im spiral arm, possibly indicating dense pre-OB star 
formation gas (see Section 2]) . 
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Fig. 2. — Locations of high and low ^2— i/i— o with respect to the molecular spiral arms and star forming regions, (a) Map of fi2-i/i-0- 
The same as Figure[l};, but with a different color palette to separate the gas with i?2_i/i_o <0.6, 0.6-0.8, and >0.8. This figure includes 
the outer areas where a large number of pixels {> 10) need to be averaged for reference, (b) Map of high i?2-l/l-0 (> 0.8). (c) Ha image, 
(d) FUV image. Contours are the same as those in Figure [Tl and are at /co(i— 0) =12, 24, 36, and 48K ■ kms^^. 



4. DISCUSSION 

GMCs arc abundant even in the interarm regions in 
M51 (Koda et al. 2009). Although our 19.7" (-780 pc) 
beam typicahy contains many GMCs, we adopt a one- 
zone approximation to obtain insights into the mean 
physical conditions of the molecular gas and its spatial 
variations. That is, we assume that the GMCs within a 
resolution element share the same average properties and 
that CO(l-O) and C0(2-l) are emitted from the same 
physical regions. 



4.1. Large Velocity Gradient (LVG) Analysis 

The Large Velocity Gradient (LVG) approximation 
is often used for the line radiation transfer of in- 
terstellar molecular lin es (jScoville fc SolomonI 119741 : 
iGoldreich fc Kwanlll974l) . This approach enables treat- 
ment of the coupled radiative transfer and molecular ex- 
citation as a local problem since the line emission in dif- 
ferent volumes of gas is doppler shifted and decoupled; 
the approach is justified by the fact that thermal line 
broadening (~ 0.2-0.3 km s~^) is much smaller than ve- 
locity dispersions within GMCs (> a few kms~^). The 
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Fig. 3. — Phase diagrams of the CO(l-O) integrated intensity 
(top panel) and -R2-1/1-0 with CO(l-O) contours (bottom). The 
azimuthal angle (phase) is defined counter-clockwise from the west. 
The plot spans over the full 360 deg, starting from 60deg to cover 
the two spiral arms without phase wraps. The counter-clockwise 
gas flow in Figure [T] is toward right in these diagrams. The gas 
with high iJ2_i/i_o appears predominantly at the downstream 
edge/side of the molecular spiral arms. 

molecular excitation conditions then depend on the local 
kinetic temperature Tkin, volume density nn^ (since low- 
J CO emission is excited collisionally) and on the local 
line opacity (i.e., CO column density A^co per unit veloc- 
ity change). Thus, (Tkin, nH2, Nco/dv) determine the 
line ratio /^2--i/i-0: or conversely, an observed i?2-i/i-o 
constrains these parameters, albeit with some degener- 
acy. We have developed our own LVG calculation code 
using the newest CO-H2 collisional cross sections from 
lYang et all (IMTl . 

The CO column density per unit velocity gradient is 
typically in the range log(A^co/'^^^ [cm^^(kms^^)^^]) = 
16.6 — 17.3 based on the observed line widths and 
assum ing a standard ISM C O/H2 abundance of 8 x 
10~^ (jSchinnerer et al.) I201C0 . For example. Galactic 
GMCs have approximately constant mean surface den- 
sities averaged over t heir entire areas (~ 170M( 7 )PC~^ ; 
iSolomon et all 119871: iRodriguez-Fernandez et"all I2006D 
and a range of velocity widths ~ 4-18 km s~^, resulting in 
the above range of Nco/dv. For the most extreme sur- 
face densities in M51's spiral arms, the mass surface den- 
sity is ~ 1000M(7)PC~^, wh ile the velocity width is 50-100 
kms~^ (jKoda et al.ll2009D . yielding a similar Nco/dv. 

Figure [7] shows the parameter ranges for Tkin and rii/^ 
which yield values of i?2-i/i-o = 0-4,0.6,0.8,0.9 and 
1.0. In each case, two curves are shown corresponding 
to log(iVco/'^^') = 17.3 and 16.6. (Note that the relative 
position of the two curves is reversed for i?2-i/i-o ^ 0.9, 



compared to that for lower values of i?2-i/i-o-) Values 
of i?2-i/i-o < 0-4 and > 1 can be obtained but are 
not shown here since they are not seen in our observa- 
tions. Galactic GMCs typically have Tkin ~ 10 K and 
~ 300 cm~ ^ (logrtg, ^ 2.5) with a s catter of about 
a fac tor of 2 (jScoville fc SandersI Il987t ISolomon et al.l 
|1987[) . yielding a range in i?2_i/i_o ~ 0.4-0.8 and a mean 
observed ratio of 0.6-0.7 from Figure[7l This range is con- 
sistent with the observed i?2-i/i-o in interarm GMCs in 
the Galaxy and the bulk parts of the Orion GMCs, i.e., 
the larger region not close ly associated wit h active high 
mass star formation (pakamoto et al.lll997l) . 

Figure [7] shows that i?2-i/i-o is a sensitive probe of 
physical conditions in the bulk molecular gas. A high 
value of i?2-i/i-o ~ 0-8 requires increases in either Tkin 
or (or both) by a factor of ~2-3, relative to the stan- 
dard CMC values of 10 K and 300 cm~^. To reproduce 
the even higher ratios ~ 0.9-1.0 within an expected range 
of Tkin or , Tkin and nH2 need to increase more than 
factors of ~3 and ~10, respectively. 

_R2-i/i-o shows spatial variations in M51. The ra- 
tio often reaches as low as ~ 0.4-0.6 in the interarm re- 
gions, while it becomes as high as ~ 0.8-1.0 in the spiral 
arms (Figure [H d andEj). The lower i?2-i/i-o (--0.4- 
0.6) generally seen in the interarm regions implies that 
the gas there has significantly lower Tkin and rifja than 
that in the arms. Based on the similarity in i?2-i/i-0: 
we conclude that the GMCs in M51's interarm regions 
have typical physical conditions similar to their Galac- 
tic counterparts. In the Galaxy, the mass dividing line 
for 50% of the total II2 content is ~ 4 x 10^ M©, i.e. 
GMCs more massive tha n this account for half th e to- 
tal molecular gas mass (|Scoville fc SandersI 119871) . In 
M51, we can derive a similar CMC mass 50% percentile 
of ~ 5 X 10 ^ MfT) using higher reso lution interferomet- 
ric imaging (|Koda et al.l 120091 [20Tl[) which is capable of 
resolving the emission of more massive clouds and com- 
paring their total emission with that in the single dish 
maps which also include contributions from lower mass 
clouds. 

The higher i?2-i/i-o '^O. 8-1.0, seen typically on the 
leading (downstream) edges of the spiral arms, requires 
higher values, both in Tkin and riua by a factor of 2-3, 
compared to those of typical GMCs and interarm GMCs. 
This downstream gas is associated with ongoing star for- 
mation (Section [S]). Therefore, increases in both Tkin 
and rifja , rather than in either Tkin or alone, may 
be more likely. On the other hand, the gas with high 
i?2-i/i-o '^0.8-1.0 on the upstream edge typically shows 
low star formation activity (Figure[S]), possibly explained 
by increases in alone, without significant increase in 

Tkin- 

Of course, this one-zone analysis is simplistic for re- 
gions extended over ^ 780 pc, but it still provides an in- 
dication that the bulk properties of the gas evolve from 
the interarm regions through the spiral arms. The inter- 
arm gas is mostly in GMCs without active star forma- 
tion, and upon entry into the spiral arms, the molecular 
gas becomes warmer and/or denser by factors of at least 
'^2 — 5, possibly due to dynamical compression (e.g., at 
the upstream edges), but primarily due to young stellar 
heating (at the downstream edges). 
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Fig. 4. — Spatial variations of i?2— i/i— o- Left: Definition of spiral arm phase (ip). A logarithmic form of spiral pattern with a pitch angle 
of 20 deg is adopted. The spiral phase starts from the thick solid curve in cyan (ip = deg) and increases counter-clockwise with thin cyan 
curves drawn at a 20 deg interval. The background image is the same CO(l-O) image as Figure [iji. Rough areas of two molecular spiral 
arms are hatched in blue {ip = 60-100 and 230-270 deg). The leading edges of the molecular spiral arms (and comparably, star-forming 
spiral arms in Ha and FUV) are hatched in red = 80-120 and 250-290 deg). A molecular spiral arm deviates from this simple logarithmic 
form at the top-right corner (indicated by gray solid lines), and we exclude this area from this analysis. The central region is defined as 
the yellow circle with a radius of 45". Right: Average -R2-1/1-0 (top) and fractions of high ratio gas (0.8-1.0; middle) and low ratio gas 
(0.4-0.6; bottom) as functions of the spiral phase. Example histograms are presented in Appendix |X] The typical standard deviation is 
~ 0.11 in J?2— 1/1— (top panel). The Poisson errors in fractions are about ~ 1% (middle and bottom). The values are calculated over 
40 deg segments in i/i. The molecular spiral arms and their leading (downstream) edges are hatched in blue and red, respectively, as in the 
left image. 



4.2. Comparison with Galactic Studies 

A similar evolution ol molecular gas, i.e., warmer 
and/or denser in spiral arms, has been observed in the 
Galaxy, although uncertainties re main due to our in - 
plane viewpoint for the Milky Way. I Sanders et al.l H1985[ ) 
found GMCs with high brightness temperature mostly 
associate d with spiral arms as delineated by radio HII 
regions. iSawada et al.l ()2012atf bl) recently demonstrated 
that the molecular gas with high brightness tempera- 
ture is located continuously al ong the spiral a r ms in 
the l ongitude-velocity diagram. ISakamoto et al.l ()1995l 
[T997I) found, from their Galactic Plane C0(2-l) survey, 
that molecular gas with the highest i?2-i/i-o is predom- 
inantly along the Sagittarius and Scutum spiral arms. 
Both the Galaxy and M51 exhibit similar variation in 



R 



2-1/1-0 



across their disks. 



The great advantage of the Galactic studies is their 
high linear reso lution. ISakamoto et al.l ()1997D and 
iHasegawal ()1997D classified molecular gas into three cat- 
egories based on their i?2- 1/1-0 measurements: low 
ratio gas (i?2-i/i-o < 0.7), high ratio gas (0.7-1.0), 
and very high ratio gas (> 1.0). The low ratio gas 
is often observed in large, n on-star forming parts of 
GMCs (ISakamoto et Il994l ) and in interarm GMCs 
(ISakamoto et al.l 11997( 1 The very high ratio gas ap- 
pears exclusively i mmediately arou nd HII regions (e.g., 
the M42 nebula; ISakamoto et aLl irggj). and hence is 
likely strongly heated by radiation at the boundary of 
the HII regions. The high ratio gas can also be found 
around the HII regions, just outside the very high ratio 
gas. In Galactic GMCs, high ratios are also associated 
with dense core regions prior to efficient star formation. 



but where the CO levels have become fully t hermalized 
and the temperatures are modestly elevated (jHasegawal 
IT997I) . 

The central ~ 2.8 kpc region of M51 shows high 
i?2-i/i-o(~ 0.8-1.0). This is similar to the Galac- 
tic center, where i?2-i/i-o ~ 1-0, and sometimes > 
1, over its central 900 pc region. These ratios are 
quite high in comparison with typical GMC values, 
though they are around the average for galactic centers 
(some starburst galaxies show i?2 -i/i-o > 1-0 their 
centers; [B raine fc CombesI [1991 ISawada et aLll200l . 
ISawada et al. (2001) explain that the high and very high 
ratios in the Galactic center are caused by intrinsically 
high Tkin and nuai combined with low opacity in low J 
transitions due to the large velocity dispersion and a very 
high excitation (J > 3). The high ratios (0.8-1.0) in the 
center of M51 probably indicate similar physical condi- 
tions to those of the Galactic Central Molecular Zone. 



4.3. Causes of High Ratios 

The spiral arm and galactic nucleus GMCs with higher 
i?2-i/i-o are likely heated primarily by associated OB 
star formation (Figure [TJ ^ and [5]), and their densi- 
ties raised by dynamical compression within the spiral 
arms (i.e., close-enco unters and/or coll isions with other 
GMCs) and shearing (jKoda et al.ll2009[ ). ah of which can 
lead to the high ratio gas (> 0.7; and up to 0.8-1.0). 

A substantial fraction of ionizing photons leaking out 
of HII regio ns jFerguson et al.l 119961: iWang et al.l 119991: 
iBlanc et al.ll2009D might be a cause of the elevated heat- 
ing outside the HII regions. However, to alter the tem- 
peratures over the very large regions (> 780 pc) at the 
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Fig. 5. — Correlations of the line ratio R2-1/1-0 with (a) the Spitzer 24fim surface brightness /24/jm, (c) CO(l-O) integrated intensity 
7qo(1— 0)1 ^tid i^) ratio of f2Aiim. and Jco(l-O)- (^1)1 (d), and (f) are the same as (a), (c), and (e), respectively, but with binned 
averages and standard deviations. The /24/jm, ^co{l— 0)5 s-^d /24/jm/Jco{l— 0) approximately proportional to the surface densities 
of star formation rate, molecular gas surface density, and star formation efficiency, respectively. Note that there are some dependences 
between x and y axes in panel (c),(d),(e), and (f), since both axes include /co{l— 0)- However, the dependences do not explain the absence 
of data in the lower-right corners, since the detected parameter ranges can cover the corners. 



leading sides of the spiral arms (as seen here), it seems 
more likely that the radiative heating is provided by the 
longer wavelength continuum of the young star clusters 
(i.e., longward of the Lyman limit) combined with hy- 
drogen recombination line radiation from the HII regions 
(since the Lyman continuum is almost surely absorbed 
over these long paths) . An abundant population of active 
HII regions exists along the spiral arms, so me of them 10- 
100 t imes brighter than the Orion nebula (jScoville et al.l 
[2OOI and of course there is considerable luminosity long- 
ward of the Lyman limit. High ratios over a large area 
are also found around 30 Doradus complex in the Large 
Magellanic Cloud, one of t he most active st ar forming 
regions in the Local Group (jSorai et al.ll2001[ ). 



We note that cloud-cloud tidal interactions within the 
spiral arms can also inject a large amount of internal en- 
ergy into GMCs in very close encounters. The typical 
gas thermal energy in a GMC is ^ 2 x 10^^ ergs; this 
is two orders of magnitude smaller than the energy in- 
crease due to a close cloud-cloud encounter, ~ 3 x lO'^^ 
erg [e stimated using the impulse approximation (jSpitzeii 
Il958f l. an en counter speed of 10 km s~ ^ , and typical GMC 
parameters (jScoville fc Sanderslll987[ )]. Dynamical envi- 
ronments in the spiral arms can possibly explain high 
ratios at the spiral arm entry (at their upstream edges). 



5. IMPLICATIONS TO OTHER STUDIES 
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Fig. 6. — Comparisons of the CO J = 2 — 1/1 — line ratio i?2-l/l-0 (Isft), Ho image (right), and contours (white) from the Spitzer 
24/im image (without smoothing) at 4, 8, and 12 MJy/str. The color scale of the left panel is the same as that in Figure [T] The high 
R2_i/i_o regions coincide mostly with star forming regions (with high /24fim and HII regions). The contours in magenta are the regions 
of high iJ2— 1/1— Oi > and low star formation efficiency, /24;tm/^co(l— 0) < ".6 (see Figure^. These regions appear both on the spiral 
arms and their trailing/upstream edges (e.g., in the north-west quadrant), possibly indicating the growth of pre-star forming dense cores. 




?.0 2.5 3.0 3.5 4.0 

losiufj^ [cm.-^]) 

Fig. 7. — Large Velocity Gradient (LVG) calculations. The line ratios -R2-i/i~o are plotted as functions of the gas kinetic temperature 
Tjjin and H2 density njjj. Most GMCs are in the range of the CO column densities log{NQQ/dv) = 16.6 (blue) and 17.3 (cyan) assuming 
the CO fractional abundance to H2 of 8 X 10~^. The regions between the curves for the two column densities (blue and cyan) are shaded 
for clarity. 
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Our results indicate that CO (2-1) emission is sensitive 
to the change of physical conditions of bulk molecular gas 
in GMCs, in contrast to CO(l-O), which traces the bulk 
molecular mass. The spatial variation of i?2-i/i-o indi- 
cates the evolution in physical conditions of the gas as 
it passes from interarm to arm regions during the course 
of azimuthal orbital motion. The use of C0(2-l) is rea- 
sonable as a tracer of molecular gas mass when galactic 
structures are u nresolved, e.g., in s tudies of high redshift 
galaxies (e.g., iFraver et al.l 120081 ) . since the -R2-1/1-0 
varies only by a factor of ~ 2. However, spatially- 
resolved studies of CO (2-1) in nearby galaxies may be 
quite misleading with respect to inferred variations in 
the gas mass surface density due to changes in -R2- 1/1-0 
with both radius (e.g. in the nuclear regions) and with 
azimuth due to the spiral arms. The systematic trends il- 
lustrated in Figure [5] will affect correlation studies of the 
gas mass and other parameters such as the Schmidt rela- 
tion (iBigiel et all 120081 : iLerov et all 120081 : iSchruba et all 
1201 It) . In fact, the Schmidt relation using the CANON 
CO(l-O) data of ^10 galaxies sh ows a non-linear, power- 
law relationship ()Liu et al.ll2011l Momose et al. in prep.), 
as opposed to the linear relation derived from the above 
quoted CO (2-1) studies. 

6. SUMMARY 

New multi-beam receivers and OTF mapping tech- 
nique in millimeter astronomy enable us to map the large 
molecular disk of M51 in CO J = 1 — and 2 — 1 transi- 
tions with unprecedented accuracy. We find systematic 
spatial variations in the line ratio i?2 -1/1-0 ~ from low 
(< 0.7; often ^ 0.4 — 0.6) in the interarm regions to high 
(> 0.7; ~ 0.8— 1.0) in the spiral arms, especially on their 
leading (downstream) edge. The central ^ 2.8 kpc region 
also exhibits high i?2-i/i-o (~ 0.8 — 1.0). 

The interarm gas is mostly in GMCs without active 
star formation, and upon entry into the spiral arms, the 
molecular gas becomes warmer and/or denser by factors 
of ~ 2 — 3. The gas in the spiral arms and galactic nucleus 
are likely heated primarily by associated OB star forma- 
tion (especially their continuum radiation longward of 
the Lyman limit) and their densities raised by dynamical 
compression within the spiral arms. Some high ratio gas 
exists on the upstream edges of the spiral arms, perhaps 
indicating gas becoming dense prior to star formation. 
Similar variations in the CO line ratios have been seen 
in the Galaxy - GMCs with inactive star formation have 
low i?2-i/i-o {"^ 0-4 — 0.6) and ones with active star 
formation show high i?2-i/i-o (^ 0.8 — 1.0). 
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AST-1211680 and NASA through grant NNX09AF40G, 
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Fig. 8. — Sample histograms of _R2— i/i— o spiral arm regions (i.e., the spiral phases of 100 deg and 260 deg in Figure |4]l, interarm 
regions (40 deg and 180 deg), the entire disk, and the central region. The average and standard deviation are in the parenthesis at the 
top-left corners. Two histograms are drawn for central regions with the diameters of 90" (~ 3.6 kpc; blue) and 70" (~ 2.8 kpc; green). 



APPENDIX 
A. HISTOGRAMS OF R2_i/i_o 

Figure [5] presents sample histograms of i?2-i/i-o in spiral arms, interarm regions, entire disk, and central region 
defined in Figure |4l Differences in the average -R2-1/1-0 between the spiral arms and interarm regions are significant 
compared to the standard deviations. 

The average ratio of the spiral arm of 260 deg is smaller than that of 100 deg. The fraction of high ratio gas is also 
lower for the 260 deg arm (Figure |4]), though the increase from interarm regions is still evident. These differences are 
partly due to the high ratio gas appearing upstream this spiral arm at smaller radii and downstream at larger radii 
(Figure [Hand [5]). A full discussion of the nature of this spatial variation (upstream vs downstream) is beyo nd the scope 
of thi s paper, but possible causes include the scattered distribution of star forming regions al ong this arm ([Egusa et al.l 
I2009D . cloud-cloud interactions, or dynamical gas compression due to the wiggle instability (|Wada fc Kodall2004[) . 

B. CORRELATIONS BETWEEN Ra-i/i-o AND STAR FORMATION RATE TRACERS 

We adopted the 24/im image as an approximate tracer of star formation rate and showed the correlation be- 
tween i?2-i/i-o star formation activities in Similar correlations exist with other tracers of star formation 
rate. Figure |9] shows cor relations with Ha flux and with star formation rate calculated by combining 24/Ltm and Ha 
(|Calzetti et"aLll2007t iKCTinicutt et al.ll200'^ . 
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